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Abstract. Quasi-hydrostatic cooling of the intracluster gas is studied. In the quasi-hydrostatic model, work done 
by gravity on the inflow gas with dP ^ 0, where P is the gas pressure, is taken into account in the thermal balance. 
The gas flows in from the outer part so as to compensate the pressure loss of the gas undergoing radiative cooling, 
but the mass flow is so moderate and smooth that the gas is considered to be quasi-hydrostatic. The temperature 
of the cooling gas decreases toward the cluster center, but, unlike cooling flows with dP = 0, approaches a constant 
temperature of ~ 1/3 the temperature of the non-cooling ambient gas. This does not mean that gravitational 
work cancels out radiative cooling, but means that the temperature of the cooling gas appears to approach a 
constant value toward the cluster center if the gas maintains the quasi-hydrostatic balance. We discuss the mass 
flow in quasi-hydrostatic cooling, and compare it with the standard isobaric cooling flow model. We also discuss 
the implication of M for the standard cooling flow model. 
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1. Introduction 

The intracluster gas, i.e. hot gas in clusters of galaxies, 
which undergoes radiative cooling is thought to flow in 
toward the cluster center to maintain pressure equilib- 
rium (see Fabian i99Hfor a review and references therein). 
According to the standard cooling flow picture, the mass 
deposition rate M can be estimated from the differen- 
tial luminosity. X-ray observations showed that M inferred 
from the surface brightness profile was over 100 Mq yr _1 
for so called cooling flow clusters. However, firm evidence 
has not been found for such a large amount of cooled gas 
that should be detected at longer wavelengths if it existed 
in the cluster core (e.g., Edge[2001). 

ASCA and recent Chandra and XMM-Newton obser- 
vations (e.g., David et al. 120011 Kaastra et al. 120011 Allen 
et al. 120011 Peterson et al. I2003f) with better spatial and 
spectral resolution reveal, however, that M is smaller by 
an order of magnitude than what was estimated before 
from the surface brightness. Moreover, spectroscopic anal- 
ysis shows that the gas temperature decreases toward the 
cluster center, as expected for cooling flows, but down to 
only half or one-third the temperature of the non-cooling 
ambient gas. This implies the lack of soft X-ray emission, 
or no appreciable contribution of lower temperature gas. 



To explain the observations, various theoretical mod- 
els have been proposed (see e.g., Bohringer et al. 2002 
Peterson et al. 120031) . Many of those consider the cen- 
tral active galaxy (with an active galactic nucleus (AGN)) 
as a heating source to compensate radiation loss (e.g., 
Churazov et al. 120011 Ruszkowski & Begelman |2I 
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Kaiser & Binney 2003). This idea is interesting from the 
aspect of history of activity (Reynolds, Heinz & Begelman 
2002}, since many cooling flow clusters possess cD galax- 
ies in their central regions. However, some fine tuning is 
required. While the radiative cooling rate is determined 
by the local values of temperature and density, the con- 
ductive or convective heating rate including hot bubbles 
from AGN depends on the scale length of the temperature 
or pressure as well. This could make the processes compli- 
cated, and the parameters such as diffusion coefficient or 
mixing length remain open (Narayan & Medvedev 120011 
Voigt et al. BlEIZt - 

In a standard analysis of X-ray data, an isobaric cool- 
ing flow model is applied to estimate M (Fabian [1994). 
However, analysis based on surface brightness and analysis 
based on spectrum give inconsistent results, as mentioned 
above, in spite of both being based on the same physi- 
cal concept. This may suggest some problem underlying 
the standard cooling flow model. In the standard isobaric 
model, work done by gravity on the inflow gas is not prop- 
erly taken into account. This is not an additional heating 
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source but is naturally expected for the gas under gravity. 
If the inflow is induced and regulated smoothly by local ra- 
diative cooling under gravity, quasi-hydrostatic structure 
may be attained. Such a simple picture is thought to be a 
starting point for considering an alternative heat source. 

In the present paper we examine quasi-hydrostatic in- 
tracluster gas undergoing radiative cooling, taking proper 
account of the work done by gravity. In Sect. 2 we present 
our quasi-hydrostatic cooling model. We compare our 
model with the standard isobaric cooling flow model in 
Sect. 3, and give concluding remarks in Sect. 4. 

2. Quasi-hydrostatic cooling model 

The internal energy U of the gas changes as 

dU = C v dT = d'Q - Pd- (1) 

P 

where Cy is the specific heat at constant volume, and T, 
p and P are the gas temperature, density and pressure, 
respectively. The first and second terms on the right-hand 
side represent heat flowing into and work done on the gas, 
respectively. If the gas is flowing toward the cluster center, 
work by gravitational force would be done on the gas. 

We consider a spherically symmetric structure. By us- 
ing the EOS (equation of state) of an ideal gas and the 
hydrostatic balance equation 

ldP=-™£dr, (2) 

Eq. Q is rewritten for the quasi-hydrostatic intracluster 

gas, as 

k fur CM* 

^-k dT = -Ldt - ^-^-Mdr. (3) 

2 /ira r z 

Here M(r) is a gas mass undergoing radiation loss at 
a rate T(r); one can consider, for instance, a shell as 
M(r) = 4irr 2 pAr. M* is the gravitational mass contained 
in the radius concerned, p the mean molecular weight, m 
the proton mass, and other quantities have their usual 
meanings. The term 'quasi-hydrostatic' does not mean 
static but is used for the steady state under gravity. 
From Eq. J2J) and the continuity equation, 



M r 



2 dr 
-Attt p—r, 
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we obtain 



5 k dT GM; 
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(5) 



where n — pj pm is the number density and n 2 A the cool- 
ing rate per unit volume. M r represents the gas mass flow- 
ing from the outer part into r per unit time. 

We assume that the intracluster gas is isothermal be- 
fore cooling when virialized, i.e., (3/2) (fcT//im) ~ —<fi 
for the gravitational potential <f> in the cluster (see e.g., 
Sarazin 1986). In such a case, the gravitational mass M*, 
which includes the galaxies and the dark matter as well 



as the gas, can be replaced with the equilibrium (virial) 
temperature To, as 



gm: 



2°* 



3 M 

r\P 1 
I /i?77 



(6) 



where o> is the line-of-sight velocity dispersion of the 
galaxies in the cluster; f3 is a parameter in the so-called 
/3-model of intracluster gas. We consider that M* is large 
enough compared to the gas mass M r within r and that 
the gas inflow hardly alters the gravitational potential. 
Eq. © with To is expressed in the form 



d\uT 
dlnr 



47r 3 n A pm 
Y r (3/2)kTW r 
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where, in the last expression, 

All n ~ 

M r = —r p and t co . 



(3/2)nfcT 
n 2 A 



(7) 



(8) 



Hereafter the quantities with tildes represent those evalu- 
ated with the local values of density and temperature; M r 
and t coo i are the mass of a uniform gas sphere and the cool- 
ing time, which are given by p(r) and T(r). Accordingly, 
M r represents the gas mass within r that cools radiatively 
in the cooling time t coo i given at r. Eq. 101 or shows the 
relation between variations dT and dr to satisfy Eq. for 
the relation between M r and L under quasi-hydrostatic 
balance, i.e., hydrostatic balance (Eq. (J2J) with a mass 
flow (Eq. Q). 

We consider that the gas inflow across r is caused by 
radiative cooling within r under quasi-hydrostatic balance 
at every r. While the gas of mass M r inside loses its pres- 
sure on a time scale (dlnP/dt)^ 1 ~ (<91nT/<9t) -1 = t coo i, 
gas of temperature T comes into the region at rate M r to 
compensate the loss with the time scale t coo i. The inflow 
mass at r during t coo i is expressed approximately as 



■ - dM r ~ 

M r t coo i ^ — tcooi rs "' JW r . 

at 



(9) 



Therefore, we assume M r ~ M r /i CO ol f° r a local quasi- 
static change at a given radius. This condition for observ- 
able quantities and its validity are discussed in detail later. 
It should be noted that M r differs from the gas mass M r 
contained within r, as it depends on the density profile 
p(r) . The cooling time of the gas M r is also different from 
t C ooi- We mention later such global quantities in the dis- 
cussion of the above condition. 

The gas that undergoes radiative cooling has a temper- 
ature profile given by Eq. J2J, while the ambient gas that 
is free from radiative cooling remains isothermal. We de- 
fine a radius r coo i at which t CO oi ~ tn — Hq 1 , the Hubble 
time, tcooi < t H at r < r coo i, and t coo i > t H at r > r coo i, 
because t coo i becomes monotonically shorter toward the 
inner region. At r > r coo i ~ 100 kpc, radiative cooling has 
little effect on the hydrostatic structure, and the gas tem- 
perature is nearly equal to the gravitational temperature. 
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On the other hand, at r < r coo i the effect of cooling be- 
comes important. The temperature decreases and the gas 
is flowing inward, yet the gas is likely quasi-hydrostatic 
if the inflow is smooth and its speed is much slower than 
the local sound speed oc VT~, or more correctly, the ram 
pressure is small compared to the thermal pressure. 

Thus, the temperature of the quasi-hydrostatic cooling 
gas begins to decrease from To at r coo i , and approaches an 
asymptotic value given by 



with the inflow velocity oc r 1 *. Here 1 — £ > in the tem- 
perature range of interest and it is likely £ < for smooth 
inflow. For a quasi-static change at given r, that is, the 
condition M r ~ M r /t coo i means 



l + tv-77(l-£)-C~fJ 



(13) 



3 M r /t c 



with a small value of |£| (= r/L v ) compared to unity, 
where L v is the scale length of the inflow velocity. 

Thus, when the gas is cooling with quasi-hydrostatic 
balancing, the density varies as 



(10) 



p oc r 



(14) 



toward the cluster center. In Fig. 1 we show a temperature roughly as p oc r ■ at r < r cool and approaches oc r 

for rj — > toward the cluster center. Note that the local 
cooling time is independent of the temperature/density 
profile, 




r/r„ 



Fig. 1. The temperature of cooling quasi- hydrostatic gas, 
as a function of r/r coo i for the case M r = f3~ 1 M r /t CO oi- 
This corresponds, for instance, to M r = M r /t CO oi for (3 = 1 
and to M r = l.5M r /t coo i for f3 = 2/3. 

profile given by Eq. Q with M r ~ M r /t CO oi at every r. For 
M r = f3~ 1 M r /t CO oi, the temperature decreases as T oc r 6 / 5 
at T - T , as T oc r 3 / 5 at T ~ (1/2)T , and eventually 
becomes ~ (1/3)T . In practice, at the cluster center, the 
gas accumulates and cools further to a lower temperature. 
It is likely that at T < T (r < r coo i) the temperature is 
somewhat less steep than T oc r 6 / 5 and becomes flatter 
outward so that the temperature connects smoothly with 
T ~ Tq at r ~ r coo i. If electron conduction works, the 
temperature would vary more gradually at 0.5 ^ r /r coo i & 
1. 

We examine the density profile of the gas in quasi- 
hydrostatic cooling considered here. For simplicity, we ex- 
press the density and temperature profiles as p oc r a and 
T oc r 71 , where rj is given by Eq. Q. Since 



M r oc r 3+a and t cool oc r^+^-V 
with A oc T*, we have 



Mr/l 



M r 



co ° l oc r 1 + Q ~ I '( 1 -?)-C 



(11) 



(12) 



^■cooi y 



i-C 



(15) 



and varies roughly as t coo i oc r. Inversely, if observations 
show t coo i oc r, this implies — a + rj(l — £) ~ 1 and therefore 
the validity of the assumption M ~ M r /t coo i at every r or 
(M r /tcooi) /M ~ constant. A recent report on 16 clusters 
observed by Chandra (Voigt & Fabian 200|J seems to sug- 
gest tcooi oc r. These properties as well as the temperature 
profile predicted by the quasi-hydrostatic cooling model 
with M r ~ M r /tcooi can be confirmed by observations. 
The gas mass M r contained within r is larger than M r for 
the density profile with a < 0. With t CO oi being the global 
cooling time of the gas M r , as 



M r 



Awr pdr 



3 + a 



-M r 



and 



_ f 4irr 2 (3/2)nkTdr 3 + 2a + ?y£ - 

tcOOl ~. i T 0~A 7 ^ ~[ | | t C Ool ) 



J 47rr 2 n 2 Ac?r 



3 + a + 77 



the difference between M r /t CO oi and M r /t CO oi is estimated 
to be within a factor of ~ 2. 

While the gas is cooling, Af r is coming in from the 
outer part so as to maintain the quasi-hydrostatic balance. 
Such a flow is not very drastic, since the gas is cooling 
monotonically and thereby the inflow is smooth. Thus, 
the mass flow rate follows 



■ 6 - ion 



dhxT 
dlnr 



3 To 
h' T 



(16) 



where L r = (47r/3)r 3 n 2 A is given by the density and tem- 
perature at r. M r is in the range (2/3-2)-L r (pm/kTo) for 
To and does not vary very much through the flow (see 
Sect. 4). For M r ~ M r /t cooh Eq. JTBJ reduces to 



Mr 



2 -L, 
3 



pm 
~kT 



8.4 



Tcool 



(lOO kpc / ' 10 ' ; cm 



(kT \ 
\keV J 



A 



10 



-23 



erg s 



Mq yr" 



(17) 



(18) 



4 



K. Masai and T. Kitayama: Quasi-hydrostatic intracluster gas under radiative cooling 



where p, = 0.61 is assumed. 

For the values of typical cooling flow clusters, the mass 
flow rate may be ~ IOOMq yr _1 . The local inflow rate M r 
does not mean the mass deposition rate (see Sect. 3) but 
merely gives an upper limit on the cooled mass in the clus- 
ter. Eq. H16| l or (|17|) means that following the change in T 
(Eqs. J3J and J5|) M r increases with increasing L r . This is 
reasonable if the mass inflow is induced by radiation loss; 
M r is not simply proportional to L r but depends also on 
T. The dependence in Eq. (|17JI should be compared with 
that in the standard isobaric cooling flow model, which is 
discussed in the following section. 



3. Comparison with isobaric cooling flow model 

In the quasi-hydrostatic model, work done by gravity on 
the intracluster gas is taken into account. The pressure is 
given by Eq. (fJJ, and therefore dP ^ in this model. In 
comparison with quasi-hydrostatic cooling, we discuss iso- 
baric cooling which is often applied to estimate the mass 
flow rate from X-ray observations. 

For given mass M, with dP = Eq. (TJJ gives 



5 M 

-k dT = d'Q = -Ldt. 

2 fj,m 



(19) 



This is nothing but the definition of the specific heat Cp — 
(d'Q/dT)p at constant pressure. Eq. (fTT)|) simply expresses 
that for mass M the temperature decreases with time due 
to radiation loss under dP = 0. The density must increase 
to maintain the pressure equilibrium, and the cooling time 
becomes shorter as cx T 2 / A. Thus, the gas cools and flows 
in increasingly toward the cluster center. 

Eq. 119|) may be compared with the expression for the 
standard cooling flow model (Fabian 1994), 

— k CF dT = dL, 
2 fim 

where dT is defined for a change from T to T — dT 
(Johnstone et al. lT^I Fabian |T3MJ> . i.e. dT > for cool- 
ing, and corresponds to the variation — dT in the present 
paper. So that, instead of the above equation, we consider 



h^dT 
2 fj,m 



-dL 



(20) 



as the standard cooling flow model. With increasing radi- 
ation loss (dL > 0), the temperature decreases (dT < 0) 
for given Mcf (> 0); otherwise, if Mcf < were the 
case, the temperature would increase for dL > 0. 

According to the standard cooling flow model (Fabian 
I1994|l . dL in Eq. IT2UI) is the differential luminosity defined 
as dL = (dL/dT)dT. Therefore, Eq. (FJUJ) is not a deriva- 
tive of Eq. ljT§|) . This implies that Eq. (T2U|) comes from 
some other concept than dP = in Eq. (Tj^) . Mc f may be 
interpreted as the gas mass that is cooling and dropping 
out of the flow per unit time (see Wise & Sarazin ri993[) 
or the mass deposition rate (Johnstone et al. I1992|) . 



For the enthalpy of the cooling gas of given T, with 
dP = we may have 

1*±T~-L, (21) 
2 [im 

consistently with Eq. (|19f) f° r given M. Eq. 121|) means 
M < (cf. Mcf > 0) or that the gas of temperature T 
decreases at a rate cx L. The lost gas must be so cool, 
much cooler than T, that it no longer contributes to L at 
all. 

Suppose that the intracluster gas has temperature T 
to emit L. With Md being the cooled mass, it increases at 
a rate roughly as Md ~ — M > in the cluster. So that, 
for the mass deposition rate M d , we have a global relation 



5 M d 

-k—T^L, 
2 fim 



(22) 



where the left-hand side does not mean the enthalpy of 
the cooling gas; remember that the gas of M d should be 
too cool (<C T) to be responsible for L. Even if Mcf is 
considered to be Md, however, it is unlikely that Eq. I|22|l 
reduces to Eq. (|20fl for thermodynamical variation dT. As 
readily seen, Eq. (|2~2"|) gives Md cx L/T, while Mcf oc 
dL/dT in Eq. I|20|l . Accordingly, the standard cooling flow 
model would yield an Mcf that traces dA/dT, which is 
determined by atomic processes regardless of the thermal 
properties of the intracluster g discussed below. 

Eq. lj2"U|) means that Mcf is proportional directly to 
—dL/dT, because dL is the variation with T in the stan- 
dard cooling flow model. Considering L cx p 2 A(M/p) cx 
pAM with P = const. = Pq, we can express the luminos- 
ity in the form 

L = ^±Po. 
\xm kT 

Thus, Eq. I|2l)fl can be rewritten as 



5 Mcf 
2 fj,m 



dL 




~dT 




U 1 - 


dlnA\ 


din I 1 / 




/ dlnA 
V 1 rflnT 


-k— 

2 fim 



(23) 
(24) 

(25) 



where 1 — d In A/d In T > in the temperature range of in- 
terest. With decreasing T, the sign of d In A/d In T changes 
to be negative at kT ~ 2 keV, below which line cooling 
dominates over bremsstrahlung. While Mcf ~ 0.5Md at 
kT £ 2 keV, M C f ~ 1.5M d at kT £ 2 keV, on account of 
the factor (1 — din A/d In T). Consequently, Mcf would 
vary sensitively to local minima/maxima in A(T). 

Eq. I|23l) has been applied to estimate the mass deposi- 
tion rate from the soft X-ray luminosity (e.g. Peterson et 
al. |5003). In a limited photon-energy range, the T depen- 
dence of A becomes complicated because of the presence of 
lines. Therefore, observations with different energy bands 
could derive different values of Mcf, or spectral analy- 
ses could derive values different from each other. In any 
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case, Mcf of the standard cooling flow model seems not 
to give a reasonable representation of the mass deposition 
rate due to radiation loss d'Q = —Ldt. If the luminos- 
ity L of a cluster is well accounted for by emission of the 
intracluster gas of nearly isothermal temperature T, Eq. 
IL'21) would give roughly the isobaric deposition rate of the 
cooled-down mass in the cluster. 

4. Discussion and remarks 

In the quasi-hydrostatic model, we consider the ideal con- 
dition that a change is quasi-static. The pressure loss 
due to radiative cooling is immediately compensated by 
continuous, smooth, spherical mass inflow from the outer 
part. Neither turbulence nor overshooting occurs during 
the process. The inflow is so moderate as not to disturb the 
hydrostatic balance significantly. Thus, if the local inflow 
rate is perfectly controlled by the local cooling rate, the 
gas is quasi-hydrostatic and M r ~ M r /t coo i is attained; 
we will investigate this point further by means of hydro- 
dynamical calculations elsewhere. 

As mentioned in Sect. 2, 'quasi-hydrostatic' means not 
static but steady-state. Accordingly, although the temper- 
ature of the cooling gas approaches a constant value to- 
ward the cluster center (Fig. 1), this does not mean that 
gravitational work can cancel out radiative cooling. This 
temperature profile is no more than a result of the steady 
inflow to maintain the hydrostatic balance. It should be 
noted that the gravitational potential or M* is assumed 
to be little affected by the gas inflow. 

At the cluster center (or close to a cD galaxy) where 
the inflow mass is accumulating, the cooling time can be 
too short for the actual flow rate to follow the cooling 
rate of the gas inside, i.e. M r /tcooi > M r . Then the quasi- 
hydrostatic balance would break, and the gas temperature 
inside decreases below ~ (l/3)/3To, yet the temperature in 
the steady flow region will continue to follow Eq. Q . 

In the opposite case, if the cooling rate at a given ra- 
dius is smaller than the inflow rate as M r /t coo i < M r , the 
flow rate M r would decrease inward and the temperature 
approaches a constant value higher than ~ (l/3)/3Tb to- 
ward the cluster center. This may be the case when cooling 
in the core is not yet very significant or when some heat- 
ing works against local radiative cooling so that t coo i is 
effectively longer than given by Eq. JSJ). 

The gas gains momentum during inflow. In the quasi- 
static model this is included into Eq. (J2J) , which determines 
the hydrostatic structure in cooperation with EOS. This is 
valid as long as the momentum flux is small compared to 
the thermal pressure. In practice, however, particularly in 
the region where the temperature approaches a constant 
value (see Fig. 1), the inflow gas may partly overshoot 
because of a fluctuation about M r ~ M r /t coo i. This may 
lead to instability and result in a break in the spheri- 
cally symmetric structure. The resultant dense part cools 
more rapidly, but the gas remains quasi-hydrostatic on the 
whole if the momentum flux is still small on average. 



Within the context of quasi-hydrostatic cooling, the 
gas inflow must be mild at every radius. However, if the 
flow evolved so that the momentum flux would become 
considerable, though still subsonic, the quasi-hydrostatic 
balance would fail unless some force works against the in- 
ertia. In such a case, for the gas to be quasi-hydrostatic, 
some momentum flux not in the form of the thermal pres- 
sure by heat may be needed against the inflow. 

Observations by coming missions can examine the ra- 
dial profiles predicted by the quasi-hydrostatic cooling 
model for several quantities, such as local cooling time 
t coo i, mass flow rate M r or local luminosity L r , and pos- 
sibly inflow velocity as well as the density and tempera- 
ture. The observed temperature profiles can be understood 
invoking quasi-hydrostatic cooling, though some momen- 
tum/heat source may be working additionally. The simple 
quasi-hydrostatic model is expected to give a basis for un- 
derstanding the thermal properties of the intracluster gas 
under gravity and for advanced study using hydrodynam- 
ics codes. 
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